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ABSTRACT 


\ 


This report summarizes the research conducted on this contract and 
previously reported on in a number of technical reports, papers presented at 
technical meetings, and published papers. The goals of the contract were to 
develop a component mode synthesis method for damped structures and to explore 
modal test methods which coula be employed to determine the relevant parameters 
required by the component mode synthesis method. 

Research was conducted on the following topics: 

1) Development of a generalized time-domain component mode synthesis 
technique for damped systems, 

2) Development of a frequency-domain component mode synthesis method for 
damped systems, and 

3) Development of a system identification algorithm applicable to general 
damped systems. 

This report presents abstracts of the major publications which have been 
previously issued on these topics. 


TIME- DOMAIN COMPONENT MODE SYNTHESIS 
FOR DAMPED SYSTEMS 

Under a previous contract, NAS8-33980, between The University of Texas 
at Austin and NASA Marshall Space Flight Center, work was begun on the 
development of component mode synthesis techniques for damped structures. 
The Final Report for that contract (1) summarizes work by Chung and Craig and 
by Howsman and Craig. The work of Howsman, which is bated on a state-vector 
component mode synthesis formulation, was continued under the present contract, 
and a report (Ref. 2) covering that research is included as an adjunct to 
this Final Report. 
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FREQUENCY-DOMAIN COMPONENT MODE SYNTHESIS 
FOR DAMPED SYSTEMS 

Time-domain component mode synthesis techniques transform the component 
differential equations of motion using component modes. Interface 
compatibility is invoked to couple the transformed component equations to 
form a coupled set of system differential equations. These can subsequently 
be solved for system modes and frequencies or integrated directly. On the 
other hand, frequency-domain analysis involves a coupling of dynamic stiffness 
matrices obtained by transforming the component differential equations of 
motion to the frequency domain by use of the Fourier transform. 

Although frequency-domain analysis is not nearly as common as time-domain 
analysis, frequency-domain techniques have been used for many years for deter- 
mining the steady-state response of structures subjected to periodic 
excitation. However, very little work has been done on the topic of determining 
the transient response of structures using frequency domain techniques. 
Since modal testing of structures frequently results in a set of frequency 
response functions describing the output/input behavior of the structures, a 
study has been conducted to determine efficient frequency-domain techniques 
for computing transient response of coupled structures. This work is reported 
on in Ref. 3, which is inc’uded as an adjunct :o this Final Report. 


SYSTEM IDENTIFICATION OF DAMPED STRUCTURES 

Under a previous contract, NAS8-33980, between The University of Texas 
at Austin and NASA Marshall Space Flight Center, Craig and Blair developed a 
generalized multiple-input, multiple-output modal parameter estimation 
algorithm (Ref. A). This frequency-domain modal parameter estimation algorithm 
has been extended by Kuidila and Craig to be applicable to systems with 
general linear viscous damping. Automated procedures for identifying 
"independent" and "dependent" measurement stations are included in the 
algorithm. The algorithm is presented in Ref. 5, which is incorporated as an 
adjunct to this Final Report. 
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ABSTRACTS OF TECHNICAL REPORTS 

A Substructure Coupling Procedure Applicable to General 
Linear Tine-Invariant Dynamic Systems (Ref. 2) 

The substructure coupling procedure presented in this report is valid 
for systems possessing general nonproportional, even nonsymraetric, damping. 
The coupled system equations of motion are derived from a variational principle, 
and f ree-interface component modes, along with a set of attachment modes 
serve as the substructure Ritz vectors, or assumed modes. A brief summary of 
the theory presented in this report follows. 

The substructure equation of motion is written in the standard matrix 
second-order form 

mx + cx + kx = f (1) 

where it is assumed that the coefficient matrices may be nonsymmetric. It 
can be demonstrated that the adjoint equation corresponding to Eq. (1) is 

my-cy+ky = f (2) 

To facilitate the solution of Eq. (1), it may be cast into state variable 
form by defining a velocity variable v through the equation 

m(x - v) ** 0 (3) 

Using variational methods, it is shown that the standard and adjoint 
differential equations may be written in state variable matrix form as 

AX + BX = F (4a) 

* 

- F 


T ’ T 
-AY + BY 


(4b) 
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where 
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One of the key motivations for substructure coupling is the possibility 
of systematicaly reducing the system order by introducing Ritz approximations 
to the substructure state variable vectors in the form 


X - $ n 
x x 


Y - 4> n (5b) 

y y 

where $ and 4> may contain, respectively, right and left substructure 
x y 

eigenvectors or other Ritz vectors. 

Coupling of substructures can be illustrated by considering two sub- 
structures a and 6 whose interface compatibility can be expressed by 


a a 
E X 


E B X B 


A similar equation holds for Y. Equations (5) and (6) may be combined, and 
the generalized coordinates separated into independent and dependent 
coordinates to produce coupling transformation matrices and such that 


(M - c T , ( nyI l 
KdI xxl IvJ 


Cn . 

y yi 


The independent system equations of motion can finally be written in the form 


An + II n « F 
s s s s s 


( 8 ) 


5 


where 


T T 

C $ A. . # C 

y ys bd xs x 

T T 

C 0 B. , # C 

y ys bd xs x 


C T $ T F 

y ys 


'xl 


and where A^j and are block-diagonal matrices 



The adjoint system equations of motion can be formed in a similar manner. 

The above substructure coupling procedure is very straightforward and is 

applicable to systems with any form of linear damping. The "price" is the 

use of state vectors and the requirement for adjoint Ritz vectors if 

biorthogonality is usea to diagonalize the substructure state matrices. 

However, the adjoint problem need not be solved otherwise. Free-interf ace 

component modes and residual attachment modes form the Ritz transformations 

$ and $ . 
x y 


Substructure Coupling in the Frequency Domain (Ref. 3) 


This report explore", procedures for coupling frequency-domain component 
models and using the synthesized frequency-domain system model to compute 
transient response. The contents of this report are briefly summarized below. 
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Let the dynamic substructure stiffness matrices be defined by 

G*(s) - m i s^+ c i s + , 1 « a,B ( 9 ) 

Therefore, the Laplace-transformed equations of motion for the ith substructure 
may be written 





( 10 ) 


where J denotes juncture (interface) coordinate, E denotes external 
coordinate, F* denotes interface forces due to adjacent substructures, and 

i i ^ 

Fj and F^ are external forces. Displacement compatibility and force 
equilibrium at the interface are given by 

a 6 

Xj - Xj - Xj (11a) 

and 

a 6 

F c + F c = 0 (11b) 


respectively. Equations (10) and (11) can be combined to form the system 
equation of motion 



r a 

g ej 


v“ 

X E 


'4 

0 4 

g ej 


X B 

X E 

• - ' 

4 

.4 4 

<4 ♦ 4 > . 


. X J . 


. r 5 * F 5 . 


Equation (12) can be written in more compact form as 


( 12 ) 


G(s) X(s) 


F(s) 


( 13 ) 
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Its solution in the Laplace domain is thus 

X(s) = H(s) F(s) ( 14 ) 

where 

H(s) - [G(s) ] -1 (15) 

Equation (14) may be employed for computing steady-state response by 
letting s » jw. However, further consideration is required if a transient 
(time domain) solution is required. In this report, Laplace-domain equations 
are considered as "offset" frequency-domain equations, and a discrete Fourier 
transform (DFT) technique is employed to compute transient solutions. 


In Eq. (13) let 


and 


X(s) ** I x(t) e~^ Ut dt 
j 0 


(16) 


F(s) - F(s) + Mxg(s) + [Ms + C] x 0 (s) 


(17) 


where 


and 


s = a + j w 


(18) 


?(s) - I f(t) e" 3t e~ Jwt dt (19) 

0 

Then, the transient response x(t) is obtained from Eqs. (.4) and (16) by 
inverse Fourier transforming Eq. (16) to obtain £(t) and then noting that 


x(t) ° x(t) e at 


( 20 ) 
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Reference 3 contains examples which illustrate the selection of 
appropriate computational parameters, such as the Nyquist frequency for 
frequency discretization and the convergence factor "a" in e -at . 


A Modal Parameter Extraction Procedure Applicable to 
Linear Time-Invariant Dynamic Systems (Ref. 5) 

In this report, a new modal parameter estimation procedure is presented. 
The technique is applicable to linear, time-invariant systems. It represents 
a further development of the multiple-input, multiple-output algorithm 
presented by Craig and Blair (4), which is an extension of Coppolino's 
Simultaneous Frequency Domain (SFD) technique (6). The principal contributions 
of the research described in this report are: 

(1) the incorporation of automatic procedures to reduce effective problem 
size , 

(2) the incorporation of general (linear) damping, rather than proportional 
damping, and 

(3) the use of more stable solution procedures than those employed by Coppolino 
(6) and by Craig and Blair (4). 

In developing a modal parameter estimation method for linear, time- 
invariant systems, it is assumed that there exists an analytical model which 
represents the dynamics of the structure and which is given by 

Mx(t) + Cx(t) + Kx(t) - f(t) (21) 

where x(t), x(t), and x(t) are the acceleration, velocity and displacement 
histories at N discrete locations on the structure. This model is reduced to 
a smaller order, which depends on the number of active modes in the frequency 
band of interest. The process of reducing, and later reconstructing, a set 
of time histories is embodiid in the following equations: 
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Y(t) 

*3 

* c x(t) 

(22) 



x(t) 

m 

\ Y(t) 

(23a) 



x(t) 

ZS 

V R Y(t) + e(t) 

(23b) 

where Y(t) 
condensing 
chapter of 

is the vector of "condensed" coordinates, 
transformation and V R the reconstructing 
Ref. 5 is devoted to automatic procedures 

, and where ¥ is the 
transformation. One 

for determining and 


Equations (21) through (23) may be combined and the resulting equation 
Fourier-transformed to give the reduced frequency domain equation 

M*Y(w) + C*Y(w) + K*Y(w) = V C D p(w) (24) 

where D is a force distribution matrix and where 

* 

M *» V MV 
C R 

C* = f c CV R (25) 

K* - » c KV r 
• • • 

Frequency-domain relationships among Y, Y, and Y may be combined wich Eq. 
(24) to give the final reduced frequency-domain equation 

Y(w) + (1/ 1w)CY - (1 /w 2 )KY - D p(w) (26) 

Least-squares solutions based on singular-value decomposition and on 
Householder transformations are employed to determine the reduced matrices C, 
K, and D. Finally, C and K are used as coefficient matrices in an eigenproblem 
whose solution gives the natural frequencies, damping factors, and mode shapes 
of the system. 
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Example problems are presented to illustrate the stability of the 
algorithm, its ability to handle non-proportional damping and ciosely-spaced- 
frequency modes, and the techniques employed for selection of the condensing 
transformation and for solution of the least-squares problem. 


ABSTRACTS OF TECBNTfAL PAPERS 

A Generalized Multiple-Input, Multiple-Output Modal 
Parameter Estimation Algorithm (Ref. 4) 

This paper, which has been accepted for publication in the AIAA Journal, 
describes a multiple-input, multiple-output frequency-domain nodal parameter 
estimation algorithm. The algorithm is restricted to systems with proportional 
damping, but is otherwise very similar to the more general procedure described 
in Ref. 5, which is abstracted above. 


Some Approaches to Substructure Coupling with Damping (Ref. 7) 

This paper, which was presented orally at the 4th International Conference 
on Applied Numerical Modeling in Tainan, Taiwan, has also been submitted for 
publication in the Journal of Astronautical Sciences. It . s condensed from 
Refs. 2 and 3, which are abstracted above, and shows the relationships between 
time-domain and frequency-domain approaches to substructure coupling. That 
relationship is summarized in Figure 1. 
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Fig. la. Time-Domain Transient Response Solution 



Fig. lb. Frequency-Domain Transient Response Solution 
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Modal Vector Estimation for Closely-Spaced-Frequency Modes (Ref. 8) 

This paper is based on a report (9) published under the preceeding 
contract, NAS8-33980. However, as part of the continuing study of Multishaker 
Modal Testing, its presentation at an international technical meeting and 
inclusion in the proceedings of that meeting are within the scope of the 
present contract. 

Complete derivations of the Asher method and the minimum coincident 
response method of modal tuning are presented in the paper. Both methods 
were applied to determine the mode shapes of a weakly-coupled structure and 
to illustrate the use of curve-fitting of experimental frequency response 
functions to obtain input data for the tuning process. The paper concludes 
that (1) the Asher method and the minimum coincident response method are 
rational procedures for employing multiple FRF's for identifying modal 
parameters, and (2) mode shapes of lightly-coupled systems are extremely 
sensitive to changes in system parameters and are therefore difficult to 
identify uniquely. 


EQUIPMENT ACQUISITION 

The following equipment was purchased for use in the research on 
multishaker modal testing covered by this contract. 

8 channels of data acquisition hardware and 

large-screen monitor to augment Genrad 2515 

computer Aided Test System $14320.15 

4 5-lb electrodynamic shakers 5736.00 

1 impact hammer including force cell 628.47 


TOTAL 


$ 20,684.62 
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CONCLUSIONS AND RECOMMENDATIONS 

Several significant contributions are contained in the reports and papers 

prepared under this contract and abstracted above. The principal contributions 

are: 

1. The development of a component mode synthesis method for systems with 
general viscous damping has been completed. By incorporating residual 
attachment modes as well as free-vibration modes, the convergence rate 
is improved. 

2. Methods of component mode synthesis in the frequency domain have been 
explored, and it has been demonstrated that digital Fourier transforms 
(DFT's) may be employed to compute transient response of a coupled 
structure. 

3. A multi-input, multi-output modal parameter estimation algorithm has 
been developed which is capable of identifying frequencies, damping 
values, and mode shapes of systems with general viscous damping. The 
algorithm is robust and requires little user interaction. 

The following recommendations are made ftr further work: 

1. Component mode synthesis procedures for damped structures which are more 
efficient than the ones currently available, including Ref. 2, should be 
developed. 

2. The application of frequency-domain techniques to compute transient 
response of coupled structures should be explored further. In particular, 
the techniques developed in Ref. 3 should be applied to systems having 
more degress of freedom. The efficiency of such a frequency-domain 
approach should be compared with time-domain methods for computing 
transient response of coupled structures. 

3. The modal parameter tstimation algorithm developed in Ref. 5 is applicable 
to general linear, time-invariant dynamic systems. It results in reduced 
system matrices, which might be very useful in an application such as 
control of flexible structures. Applications of the algorithm should be 


rw** *• 
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explored, and it should be compared with other modal paraneter estimation 
algorithms, such as the Polyreference algorithm. 
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